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Abstract On the ground of the Landau criterion we study the behavior of critical 
velocities in a superfluid two-component Bose gas. It is found that under motion 
of the components with different velocities the velocity of each component should 
not be lower than a minimum phase velocity of elementary excitations (s_). The 
Landau criterion yields a relation between the critical velocities of the components 
(v c i , v C 2). The velocity of one or even both components may exceed s_ . The max- 
imum value of the critical velocity of a given component can be reached when the 
other component does not move. The approach is generalized for a two-component 
condensate confined in a cylindrical harmonic potential. 
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^. . 1 Introduction 

q 

r — . At present considerable attention is given to the study of two-component Bose- 

Einstein condensates. The progress in cooling and trapping of atomic rarefied 
gases allows to obtain such condensates experimentally 'i 2 i 3 A 5 . From the funda- 
mental point of view, they are considered as objects in which some cosmological 
and astrophysical processes can be modelled 6 ' 7,8 . In two component systems the 
superfluid components may flow with different velocities. In such a situation an 
unusual (nondissipative) kind of the drag effect takes place 9 . It makes possible 
to create a controlled phase difference between two Bose-condensates placed in 
a two-well potential, and to observe effects 10 , similar to ones that occur in su- 
perconductive systems with Josephson contacts in magnetic fields. In this article 
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we study a related problem, namely, the critical velocities in a two component 
superfluid system with components flowing with different velocities. 



2 The Landau criterium for the two-component system 

According to the Landau criterion, the critical velocity in a single-component con- 
densate is determined by the expression 
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where Eo(k) is the excitation spectrum in an immovable condensate and k is the 
wave number. The velocity is given in the frame of reference connected with walls 
or obstacles. Eq. (fl} can be applied to a two-component condensate only in the 
case when both components move with the same velocity. 

The Landau criterium can be reformulated as the requirement of positivity of 
energies of elementary excitations in the frame of reference, connected with walls 
or obstacles. These energies depend on the velocities of the components Vi and 
\2 (in the same frame of reference), and the Landau criterium should yield some 
joint condition on Vi and v 2 . 

To find the energy of elementary excitations we use the Gross-Pitaevskii equa- 
tion for the two-component system 
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where 1//7 are the wave functions of the components, mi are the masses of the 
particles, yi, y\% are the interaction constants (yj = A%h ajj/m;, ju = 2nh (m\ + 
m2)ayij (m\m2) where a^ are the scattering lengths). 

The wave function of the component can be represented as the sum of the 
stationary part and the fluctuating part \j/j — xj/qj + 8\jfj, where 8\jfj <C 1/%. The 
stationary part of the condensate wave function can be presented in the form: 
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where jiiy 



mjyj 



' Yj n j + Y12 «3-; are the chemical potentials of the components. 



The gradients of the phases <j>y are connected with the superfluid velocities by the 
relation V; = — V<p,-. 
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The fluctuating part can be written as 
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The substitution of (|4]l into the linearized version of (f2| leads to a system of equa- 
tions for u — v coefficients, whose determinant gives the dispersion equation for 
the spectrum of elementary excitations 
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where Ej = yfF^~+2yj]nJ) is the Bogolyubov spectrum for the j'-component 
(in the absence of interaction between the components), and £ ; = h 2 k 2 /2irij. 
For Vi = v 2 = v the equation © yields 
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In this study we assume that the condition of stability of a two-component con- 
densate relative to phase separation is fulfilled (yi y 2 ~ Yn > 0). 

The requirement of positivity of (O at all k is equivalent to the condition (fl}. 
This condition yields the following expression for the critical velocity 
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where Sj = y^YjUj/nij are the bare velocities of the sound modes for the compo- 
nents (in the absence of interaction between the components). 

In a general case Vi 7^ V2 the Landau criterium requires the existence of two 
positive solutions of Eq. I© at all k. This requirement is equivalent to the following 
two inequalities: 

[E\ - (fivi • k) 2 ] [El - {hx 2 • k) 2 ] - 4e 1 e 2 rf 2 n 1 n 2 > 0, (8) 

E\ > {h\i ■ k) 2 (or El > {h\ 2 • k) 2 ) . (9) 

The critical values of Vi and v 2 correspond to the case, when the inequality ([8| 
turnes into the equality at least for one k. In the case considered the sufficient 
condition for fulfilling the inequalities ®, (O for all k is their fulfillment at k — > 
for all directions of k, complanar to Vi and v 2 . Therefore, the inequalities ©, © 
can be replaced with the system of inequalities 

(s 2 -v 2 cos 2 a) (^-v^cos 2 (0-a))-^-i 2 ^>O, i 2 >v 2 cos 2 a (10) 

(where 6 is an angle between Vi and v 2 , and a is an angle between k and Vi), that 
should be fulfilled for all a. The results of the analysis of dlOl) is given in Fig.l 

If the components move in the same direction their critical velocities are re- 
lated by the equation 

(A-vl i ){s 2 2 -vl 2 ) = ^-s\sl (11) 
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with the additional condition v c \ < s\. According to Jilt only for Vi = v 2 the 
critical velocity coincides with the velocity of the lowest hydrodynamic mode s~. 
In a general case one of the velocities may exceed i_ (the second velocity should 
be less than i_). If one of the components does not move the velocity of the other 
component may reach the maximum value 
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Fig. 1 Connection between the critical velocities at different 9 (for the parameters Yn 
v/7i 72/2). The shown value S- corresponds to the case si = S2- 



One can see from Fig.l that at d 7^ 0, n both components may move with ve- 
locities that exceed s_. Under the motion of the components in mutually perpen- 
dicular directions (6 = Tt/2) the velocities can simultaneously reach the maximum 
critical values (112b . 

Here we do not consider the possibility of excitation of vortices. Therefore, 
strictly speaking, our analysis yields only the upper bound for the critical veloci- 
ties. Nevertheless, in a number of situations the estimation for the critical veloci- 
ties presented in this paper is justified completely. For example, this occurs when 
the superfluid flows past an obstacle with a small (less than the healing length) 
linear size 11 ' 12 ' 13 . 



3 Critical velocities in the two-component condensate, confined in a 
cylindrical harmonic potential 

In Bose gases confined in optical or magnetic traps surface excitations have the 
minimum phase velocity 14 . In such systems the process of generating of vortices is 
connected with the excitation of surface modes, and the critical velocity coincides 
with the phase velocity of the lowest surface mode 15 . 

Let us study the critical velocities for a two-component Bose gas, confined 
in a harmonic cylindrical potential V(r) = mft^ 2 (i 2 +y 2 )/2. We will consider the 
case when superfluid flows are directed along z. For simplicity we assume n\ (r) — 
ni(r) = n(r) (r - radial coordinate), m\ = mi =m,y\=yi= 7 and < 712 < 7 Let 
us consider the system whose Thomas- Fermi radius 7?tf = [2(7+ Yn) n o/ m (Oo] l ^ 2 
is much large than the oscillator length of the trap. To find the spectrum of elemen- 
tary excitations we pass from the Gross -Pitaevskii equation to the linearized sys- 
tem of hydrodynamic equations for the densities nj(r,t) and the velocities Vj(r,t) 



of the components: 

- 5 ££+V(»5vy+v< v 5ny) = 0, 

L (13) 

m J + V(ySrij + Yn8ni-j + m\ 0j ■ 5vy) = 0, 

where n = «o(l — 7' 2 /^tf) i s the equilibrium density of the components (no is the 
density in the center of trap), Vqj — (0, 0, vj) are their superfiuid velocities and 8rij 
and 8\ j are the fluctuations of these values. 

The modes we are interested in are localized near the surface and the problem 
considered can be reduced to the problem for the spectrum of excitations in a Bose 
gas in a linear potentia l iW 8 . With this simplification we obtain the following 
dispersion equation 
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The region of the applicability of ( 114b is bounded from above by the condition k < 
k m where k m is the wave vector for which the contributions of kinetic and potential 
energies to the excitation energy become comparable. The value of k m can be 
estimated by equating the kinetic energy of the particles ft < k 2 /2m to the energy 
of the lowest mode (in the hydrodynamic limit) Eq- = h[2(y — yi2)nok/mRjp] 1 ' 2 

/ itinc] (y — Y\ ? ) 
what gives k m = 2? '= . At the point k w k m the dependence of the 
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excitation energy on k has a bend and the critical velocity can be estimated by 
substituting of k = k m into the equation ( 1141 ). 

As it follows from (114b for vi = V2 the critical velocity is the minimum phase 
velocity of the lowest surface mode 

l , c=S!f= (^y ,2 ,ra (15) 
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If only one component moves, the critical velocity reaches the value 
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Thus, the highest velocity can be reached when only one component flows. 

4 Conclusions 

In conclusion, let us discuss a number of possibilities to observe the predicted be- 
havior. One of them is to create barriers not penetrable for one of the components. 
If in the absence of the barriers the components move with equal velocities the ap- 
pearance of such barriers may increase of the critical velocity. The motion of the 



components with different velocities (and in different directions) can also be real- 
ized in systems, in which the components are separated spatially, for example, in a 
bilayer geometry (in that case a long-range interaction between the components is 
required for observing the effect predicted). The effect is also may be observed in 
multilayer condensates of electron-hole pairs which can appear in semiconducting 
heterostructures with an even (> 2) number of two-dimensional electron layers^. 
In electron-hole condensates the critical velocities can be easily measured since 
they are proportional to the value of critical currents. 
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